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Abstract. We introduce in this paper embedded Gaussian unitary ensemble of random matrices, 
for m fermions in Q number of single particle orbits, generated by random two-body interactions 
that are 5(7(4) scalar, called EGUE(2)-S'(7(4). Here the S(7(4) algebra corresponds to Wigner's 
supermultiplet SU (4) symmetry in nuclei. Formulation based on Wigner-Racah algebra of the em- 
bedding algebra U (451) D U (Q.) ® SU (4) allows for analytical treatment of this ensemble and using 
this analytical formulas are derived for the covariances in energy centroids and spectral variances. 
It is found that these covariances increase in magnitude as we go from EGUE(2) to EGUE(2)-s to 
EGUE(2)-5'[/ (4) implying that symmetries may be responsible for chaos in finite interacting quan- 
tum systems. 
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1. Introduction 

Hamiltonians for finite quantum systems such as nuclei, atoms, quantum dots, small metal- 
lic grains, interacting spin systems modeling quantum computing core, Bose condensates 
and so on consist of interactions of low body rank and therefore embedded Gaussian en- 
sembles (EGE) of random matrices generated by random interactions, first introduced in 
1970 in the context of nuclear shell model and explored to some extent in the 70's and 80's, 
are appropriate for these systems. Note that EGE's that correspond to the classical ensem- 
bles GOE, GUE and GSE are EGOE, EGUE and EGSE respectively. With the interest 
in many body chaos, EGE's received new emphasis beginning from 1996 and since then 
a wide variety of EGE's have been introduced in literature, both for fermion and boson 
systems [1^]. See [1,2,5,6] and references therein for recent applications of EGE's. 

EGE's generated by two body interactions [EGE(2)] for spinless fermion systems are the 
simplest of these ensembles. For m fermions in N single particle (sp) states, the embed- 
ding algebra is SU{N). It is well established that SU{N) Wigner-Racah algebra solves 
EGUE(2) and also the more general EGUE(/c) as well as EGOE(A;) [7,8]. Realistic systems 
carry good quantum numbers (for example spin S for quantum dots, angular momentum 
J for nuclei) in addition to particle number m, therefore EGE's with good symmetries 
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should be studied. EGUE(2)-s and EGOE(2)-s, for fermions with spin s = ^ degree of 
freedom, are the simplest non-trivial EGE's with immediate physical applications. For m 
fermions occupying number of orbits with total spin S a good quantum number, the 
embedding algebra for EGUE(2)-s and also for EGOE(2)-s is U{2^) D U{n) SU{2) 
[5,9]. In particular the EGOE(2)-s with its extension including mean-field one body part 
has been extensively used in the study of quantum dots, small metallic grains and atomic 
nuclei [6,10,11]. 

Wigner introduced in 1937 [12] the spin-isospin SU (4) supermultiplet scheme for nu- 
clei. There is good evidence for the goodness of this symmetry in some parts of the periodic 
table [13] and also more recently there is new interest in SU (4) symmetry for heavy N ~ 
Z nuclei [14]. Therefore it is clearly of importance to define and study embedded Gaus- 
sian unitary ensemble of random matrices generated by random two-body interactions with 
SU{4:) symmetry, hereafter called EGUE(2)-S'?7(4). Given m fermions (nucleons) in fl 
number of sp levels with spin and isospin degrees of freedom, for SU (4) scalar Hamil- 
tonians, the symmetry algebra is U{4Q) D U{Q) €5 S'C/(4) and all the states within an 
SU{'i) irrep will be degenerate in energy. Our purpose in this paper is to define EGUE(2)- 
5'C/(4), develop analytical formulation for solving the ensemble and report the first results 
for lower order cross correlations generated by this ensemble. Now we will give a preview. 

Section 2 gives a brief discussion of SU (4) algebra. EGUE(2)-S'C/(4) ensemble is de- 
fined in Section 3. Also given here is the mathematical formulation based on Wigner-Racah 
algebra of the embedding U{4Q) D U (8) SU (4) algebra for solving the ensemble. In 
Section 4, analytical formulas for m fermion U{^1) irreps fm = {4'', p} are given for the 
covariances in energy centroids and spectral variances generated by this ensemble. Section 
5 gives discussion of some numerical results. Finally Section 6 gives summary and future 
outlook. 



2. Preliminaries of ?7(4f2) D U{fl) O 5J7(4) algebra 

Let us begin with m nucleons distributed in V, number of orbits each with spin (s = i) 
and isospin (t = i) degrees of freedom. Then the total number of sp states is N = 40, 
and the spectrum generating algebra is U (40) . The sp states in uncoupled representation 
are a| ^ |0) = \i, a) with i = 1, 2, . . . , denoting the spatial orbits and a = 1, 2, 3, 4 are 
the four spin-isospin states |ms ,mt) = —^), \~^^ 5) and | — — respec- 

tively. The (4f2)^ number of operators Cia-jp generate [/(4f2) algebra. For m fermions, 
all states belong to the J7(4f2) irrep {!"'}. In uncoupled notation, Cia-jfj = o\a'^j,i3- 
Similarly U{il) and C/(4) algebras are generated by Aij and Bap respectively, where 
Aij = Yltt^i (^ia;ja and Ba/3 = X^^i Cia-ii3. The number operator n, the spin operator 
S = Sj^, the isospin operator T = and the Gamow-TeUer operator ar = {crr)^ ^, of 
U{4:) in spin-isospin coupled notation are [15], 





(1) 
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Note that aj;^,,^^ = (-l)^+''"+'"aj;_^^,_^j. These 16 operators form U{A) algebra. 
Dropping the number operator, we have SU (4) algebra. 

For the U'(4) algebra, the irreps are characterized by the partitions {F} = 
{Fi,F2,F3,F4} with -Fi > F2 > F3 > F4 > and m = Y^Li^^- Note that F„ 
are the eigenvalues of B^a- Due to the antisymmetry constraint on the total wavefunction, 
the orbital space U {Cl) irreps {/} are given by {/} = {F} which is obtained by changing 
rows to columns in {F}. It is important to note that, due to this symmetry constraint, for 
the irrep {F} each Fj < where j = 1, 2, 3, 4 and for the irrep {/} each fi < 4 with 
i = 1,2, . . . The irreps for the SU{4) group are characterized by three rowed Young 
shapes {F'} = {Fi , F^, F^} = {Fi- F4,F2- F^, F3- F4,}. Also they can be mapped to 
SO{6) irreps [Pi,P2, P3] as the 5J7(4) and 50(6) algebras are isomorphic to each other, 
lPi,P2,P3] = [(i^i+F2-F3-F4)/2, {F^-F2+F3-Fi)/2, {F, - F^ - F; + F^)/2]. 
Before proceeding further, let us examine the quadratic Casimir invariants of U (fi), C/(4), 
5f/(4) and 5*0(6) algebras. For example, 

C2 [C/(4)] = Ba,0B0,a ^ C2 [Um + C2 [U{4)] = n (f2 + 4) . (2) 

Also, in terms of spin, isospin and Gamow-Teller operators, C2 [5f (4)] = C2 [50(6)] = 
^2 _|_ y2 _|_ ^p.^^ . (^(jr). Now we have the general results, 

(O2 [t/(4)])^^> =Y.F,{F, + 5 - 2*) = /02 [5C/(4)] + ^ 

(O2 [50(6)])[^1 = (O2 [5f/(4)])^^'> = Pi (Pi + 4) + P2(P2 + 2) + Pi . (3) 

In order to understand the significance of SU (4) symmetry, let us consider the space ex- 
change or the Majorana operator M that exchanges the spatial coordinates of the particles 
and leaves the spin-isospin quantum numbers unchanged, 

M\i,a,a';j,(3,(3') = \j,a,a';i,(3,(3') , (4) 

where a, [3 are labels for spin and a',/?' are labels for isospin. As \i,a,a';j,(i,P') = 
ala,a'"'l,/3,0' Eqs. (4), (2) and (3) in that order will give, 

2M= {4,a,a4,0jr) {"la.a'aip,/}')^ 

= C2 [U{n)] - On = 4n - 02 [U{4)] 
^QM = a|2n(^l + ^) -i02[5[/(4)]| . (5) 

The preferred U{fl) irrep for the ground state of a m nucleon system is the most symmetric 
one. Therefore (O2 [?7(0)]) should be maximum for the ground state irrep. This implies, 
as seen from Eq. (5), the strength a of M must be negative. As a consequence, as follows 
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Figure 1. Schematic representation of the Young tableaux fm = {4'',p} with p = 0, 1, 
and 2. (1) = {4'-}, h = {2}, /™_2 = {4'-i, 2} and here a = 6; (2) = {4''}, /a = 
{12}, /„_2 = {4'-2, 32} and Tat = 1; (3) = {4^ 1}, /2 = {2}, /„_2 = {4'^-!, 2, 1} 
and here a = &; (4) /,„ - {4^ 1}, /2 = {2}, {P}, /„_2 = {4'-i,3} and Tab = 4; 
(5) = {4^ 1}, /2 - {12}, /„_2 = {4'-2, 32, 1} and t^c = -1; (6) U = {4^ 2}, 
/2 = {2}, f,n-2 = {4^-1, 22} and here a = 6; (7) /„, = {4^ 2}, /2 = {2}, /„_2 = {4'"} 
and here a = 6; (8) /„ = {4^ 2}, /2 = {2},{l2}, /„_2 = {4--\3,l}, Tab = 3; (9) 



= {4^ 2}, /2 = {12}, fm-2 = {4''-2, 32, 2} and Ta 



-1. 



from the last equality in Eq. (5), the ground states are labeled by SU{4) irreps with smallest 
eigenvalue for the quadratic Casimir invariant consistent with a given (m, T^), T = \Tz\. 

For even-even nuclei, the ground state SO{6) irreps are T = |N— Z|/2. For odd-odd 
N=Z nuclei, the ground state is [1] and for N 7^ Z nuclei, it is [T, 1]. For odd-A nuclei, the 
irreps are [T, i, ±i] . Therefore, for N=Z even-even, N=Z odd-odd and N=Z±1 odd-A 
nuclei the U{fl) irreps for the ground states are {4''}, {4'', 2}, {4*^, 1} and {4'', 3} with 
spin-isospin structure being (0,0), (1,0) © (0, 1), (^, ^), and (^, i) respectively. For 
simphcity, in this paper we will present final results only for these U{Q) irreps. Other 
irreps will be considered elsewhere. Now we will define EGUE(2)-S'f/ (4) ensemble and 
derive some of its properties. 



3. EGUE(2)-5J7(4) ensemble: Definition and U{4:n) D U{n) (g) 5J7(4) Wigner-Racah 
algebra for covariances 

Here we follow closely the approach used for EGUE(2)-s recently [9]. Let us begin with 
normaUzed two-particle states I/2-F2; ^'2/32 ) where the U{4) irreps i^2 = {1^} and {2} and 
the corresponding U{d) irreps /2 are {2} (symmetric) and {l2} (antisymmetric) respec- 
tively. Similarly V2 are additional quantum numbers that belong to /2 and /32 belong to 
F2. As /2 uniquely defines F2, from now on we wiU drop F2 unless it is exphcitly needed 
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and also we will use the /2 ^ F2 equivalence whenever needed. With A^(/2?;2/?2) and 
^(/2^^2/32) denoting creation and annihilation operators for the normalized two particle 
states, a general two-body Hamiltonian H preserving SU (4) symmetry can be written as 

^= E V^^^.^j{2)A^{hvfp2)A{hvlP,). (6) 

InEq. (6), V},„.„/(2) = (/2i'|/32 | | /2V2/32) independent of the /32's. ForEGUE(2)- 
5C/(4) the Vf^yiyf (2)'s are independent Gaussian variables with zero center and variance 
given by (with bar representing ensemble average), 

Thus V{2) is a direct sum of QUE matrices for F2 = {2} and F2 = {1^} with variances 
for the diagonal matrix elements and {Xf^)'^ /2 for the real and imaginary parts of 
the off-diagonal matrix elements. As discussed before for EGUE(k) [7] and EGUE(2)-s 
[9], tensorial decomposition of H with respect to the embedding algebra U{fl) SU (4) 
plays a crucial role in generating analytical results; as in |7|, the U (il) ^ SU (il) corre- 
spondence is used throughout and therefore we use U(fl) and SU (fl) interchangeably. As 
H preserves St/ (4), it is a scalar in the SU{4:) space. However with respect to SU the 
tensorial characters, in Young tableaux notation, for = {2} are F^, = {0}, {21^~^} and 
{42^~2} with u = 0,1 and 2 respectively. Similarly for /2 = {P} they are_F,^ = {0}, 
|2iii-2T. jjjjjj {221^^-4} with V ^0,1,2 respectively. Note that F^ = f^y^J^ where % is 
the irrep conjugate to /2 and the x denotes Kronecker product. Then we can define unitary 
tensors B's that are scalars in SU{A) space, 

B{f2F^0J.) = J2 ^Hf 2VIP2) A{f2vi(32) {f2vi h4 \ F.W.) 

V'2,vf,p2 

X {F2P2 FhP^ I 00) . (8) 

In Eq. (8), (/2 ) are SU{n) Wigner coefficients and {F2 ) are 6'f/(4) Wigner 

coefficients. The expansion of H in terms of B's is, 

H= Yl W{f2F^oj^)B{f2F^u;^) . (9) 

f2,Fu, UJu 

The expansion coefficients VK's follow from the orthogonality of the tensors _B's with 
respect to the traces over fixed /2 spaces. Then we have the most important relation needed 
for all the results given ahead, 

W{hF,uj,)W{f^F',ujl) = Sf2f^6p^p-J^^^- (A/j'd(F2) . (10) 

This is derived starting with Eq. (7) and substituting, in two particle matrix elements V, 
for H the expansion given by Eq. (9). Also used are the sum rules for Wigner coefficients 
appearing in Eq. (8). 
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a= 10 




(a) m (b) ™ m' 

Figure 2. (a) Widths a{m, /„) = [ (i/s)™'/- ]i/2 for = 6 and = 10 examples, 
(b) Cross correlations for = 10 examples. Note that fm = {^^,p}', m = + p and 
fm' = {4*, g}; m' = 4s + q. See text for details. 

Table 1. Pf-jm, fm) for fm = {4'',p}; p = 0, 1, 2 and 3 and {/2} = {2}, {l^}. 
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h = {2} 


/2 = {1^ 




-3r(r + 1) 


— 5r(r — 1) 


{4^,1} 


- — (2r + 3) 


-|(2-1) 


{4'', 2} 


-(3r2 + 6r + l) 


-5r2 


{4^,3} 


-5(r + 2)(2r + l) 


-|(2r + l) 



Turning to m particle iJ matrix elements, first we denote the U (O) and U (4) irreps by 
fm and respectively. Correlations generated by EGUE(2)-S'J7(4) between states with 
{m, fm) and (m' , fm' ) follow from the covariance between the m-particle matrix elements 
of H. Now using Eqs. (9) and (10) along with the Wigner-Eckart theorem applied using 
SU {Cl) (g) SU{4:) Wigner-Racah algebra (see for example [16]) will give 



= {fmFmvLP I H I fmFmVinP) {fm'Fm'vi,P' \ H \ fm'Fm'V^^P') 
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^ X] ^dtf^) ^ -^(-^^-F^) III /m>p ifm' III B{f2F^) III 

fm'V^i ) ; 

III Vt(f IT ^ III f \ - ■^f^-^ U{fmf2fm.f2; /m-2-F,y)p 

(/mill B{}2F^) III - 2^ ^ (m) — ^ — — . (11) 

Here the summation in the last equaHty is over the muhipHcity index p and this arises as 
/to (8 -Fj/ gives in general more than once the irrep/m. In Eq. (11), i^(m) — — m(m— 1)/2, 
d{fm) is dimension with respect to U{Q.) and 7V/^ is dimension with respect to the Sm 
group; formulas for these dimensions are given in [17]. Similarly, (. . .) and U{. . .) are 
SU{Q) Wigner and Racah coefficients respectively. 



4. Lower order cross correlations in EG\JE(2)-SU (4) 

Lower order cross correlations between states with different {m,fm) are given by the 
normalized bivariate moments E^r (m, fm '■ m', fm'), r = 1, 2 of the two-point function 
where, with p™'f"^ (£) defining fixed-(m, fm) density of states, 

Smfm-.m'f^, ^ pm,U(^E)p^' 'fm' {E') - p^'fm{E) p^' 'fm' [E') ; 

Sn (m, fm : m', fm') = (F)™'^'" {H)"'''^-' /^/ {H^)"^'f- ^H^)rn'j^, ^ 



S22 (m, fm : m', fm') = (iJ^)™./-^ (if2)™'./™'/ (^2)™./™ (ijzp' J. 



-1. (12) 



In Eq. (12), (^lJ'2'^"^'fm (jjg second moment (or variance) of p'^'fm{E) and its centroid 

(i?)"'^- = by definition. We begin with (H)"^'^'- (7?)"''^"' . As (iJ)™-^™ is the trace 
of H (divided by dimensionality) in {m, fm) space, only F,^ = {0} will generate this. 
Then trivially, 

(F)™'^- {H)^'^fm' = V P^^ (m, /to) (m', /to') . (13) 



h 



Note that P^^ (m, /,„) = F{m) Y.f^_^ Wf^-2/-^fJ- The formulas for pf'^ (to, /„) are 
given in Table 1. Writing ^ij2|)™^/"> explicitly in terms of to particle _ff matrix elements, 
^^2^./™ = [^(/to)]"^^,,^,,;^ and applying Eq. (11) and the 

orthonormal properties of the SU (fl) Wigner coefficients lead to 

^^^^Etttt E Q^ih-'^Jm). (14) 

/2 ^ !^=0,1,2 
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Table 2. (g2p./m^ Q'^=^'^(/2 : m, /^) and By=^{m, fm) for some examples 



fr, 



{4^,1} 
{4^,2} 

{4^,3} 



r-(n - r + 4) 



A^2}3('- + - r + 3) + A2^2j5(r- - - r + 5) 



[A|2}{6r(n - r + 1) + 90 + 15} 

+A2^2}5{2r(0 - r + 5) - n - 9} 

-^{2} 5 [S''" - + 2)r3 + (3Q2 + en - 5)r2 
+ (0 + 2)(6n + 17)r + mp. + 1)] 
+A2^,j|l(n - r + 4){(0 + 4)r - - 3} 

i [A^2}3('' + 2)(f2 -r + 2)(2rQ - 2r2 + 6r + Q + 1) 
+A2^2 j5r(n - r + 4)(2rn - 2r2 + 6r + H - 1) 



i,rx i^x 1 9r-(r + l)2(a-r)(Q + l)(n + 4) 

^ ^ 2(n + 2) 

3rn(r + l)(n -r+ l){n - r){n + i){Cl + 5) 



2 

{12} 1 
2 



4(n + 2) 

2^ ^ 25r(r - lf{n - r){n - + 4) 
2(n - 2) 

5rn(r - l)(a + 3)(n + 4)(f2 - r)(a - r - 1) 
4(Q - 2) ~ 



{4-^} 



15r /n^ 
2 V ^2 - 4 



The functions Q^{f2 '■ rn, fm) involve SU{fl) Racah coefficients and they are available 
in various tables in a complex form involving functions of Tab [18]. Here Tab are the 
axial distances for a given Young tableaux (see Fig. 1 for examples). Evaluating all the 

functions, we have derived analytical formulas for Q'^ (/2 : m, /,„) and also for ^772^™ J™ _ 
Some of these results are given in Table 2. It is easily seen that Q''"°(/2 : m, fm) = 
[pf^ {m,fm)] . Results in Tables 1 and 2 will give formulas for the covariances Sn in 
energy centroids. Similarly, analytical results for covariances S22 in spectral variances are 
derived using Eqs. (11) and (12) and then, 
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'^22{m,fm;m',fm') = 



X|2} + ^{1^} + 4X{i2}{2} 



= E [rf(^^)]"'Q''(/2:m,/™)Q''(/2:m',/„0, 

K/2)] ,.=0,1,2 

^^^^>^'> = rf({2}jd({in) ^ I^^^'^)]"' • (1^) 

Here d{F^) are dimension of the irrep F^,, and we have (i({0}) — 1, (i({2, 1^^^^}) = 
- 1, d({4,2"-2}) = 02(Q + 3)(f7 _ i)/4^ and fi({22, 1^^-*}) = n'^{n - 3)(f7 + 
l)/4. Again the functions i?" involve 5?7(0) Racah coefficients and R^^^{m, fm) = 
P{2}(^, /„)p{i'}(m, Formulas for R-=\m, /„) for = {4^p}, = - 3 are 
derived and the result for {4''} is given in Table 2 as an example. Complete tabulations for 
Qu=i,2 ajjj R'^=^ will be reported elsewhere. Equations (13), (14) and (15) are similar in 
structure to the corresponding equations for EGUE(2)-s [9]. However the functions P's, 
Q's and Rs are more compUcated for EGUE(2)-5f/(4). 



5. Results and discussion 



Numerical calculations are carried out for (iJ2j,™./™^ and 1122 for some U, — Q 
[(2slrf)-shell] and f2 = 10 [(2pl/)-shell] examples. Here we have employed A^^aj = 

A|2} = 1. Fig. 2a shows the variation in the spectral widths a{m, fm) = [ (_ff2^™'/™ ji/2 

with particle number m. Notice the peaks at m = 4r; r = 2,3, Except for this 

structure, there are no other differences between {4''} and {4'', 2} systems i.e. for ground 
states of even-even and odd-odd N=Z nuclei. Results for the cross correlations Sn and 
S22 are shown in Fig. 2b. It is seen that [En]i/2 and [1:22]^/^ increases almost linearly 
with m. At m = 4r, r = 2, 3, . . . there is a slight dip in [Sn]^/^ as well as in [S22]^/^. 
For = 6 with m = to', pii]^/^ ~ 10 - 28% and [1:22]^/^ 6 - 16% as m changes 
from 4 to 12. Similarly for m ^ to', [Sii]!/^ ~ 10 - 24% and [£22]^/^ ~ 6 - 12%. 
The values are somewhat smaller for O = 10 (see Fig. 2b) which is in agreement with 
the results obtained for EGOE(2) for spinless fermions and EGOE(2)-s. For further un- 
derstanding we compare, for fixed N, these covariances with those for EGUE(2) and 
EGUE(2)-s. Using the analytical formulas given in [7] for EGUE(2), [9] for EGUE(2)- 
s and the present paper for EGUE(2)-S'i7(4), it is found that the magnitude of the co- 
variances in energy centroids and spectral variances increases with increasing symmetry. 
For example, with = 24 [so that = 12 for EGUE(2)-s and 17 = 6 for EGUE(2)- 
SU{4)] the results are as follows. For to = to' = 6 (to = to' = 8) we have: (i) 
[Eii]i/2 = 0.017(0.026) and [^22]^/^ = 0.006(0.006) forEGUE(2); (ii) forEGUE(2)- 
s with S* = 5' = 0, = 0.043(0.066) and [£22]^/^ 0.017(0.021); (iii) for 

EGUE(2)-S';7(4), [Snj^/^ = 0.124(0.16) and [1:22]^/^ = 0.069(0.082). As fluctuations 
are growing with increasing symmetry, it is plausible to conclude that symmetries play a 
significant role in generating chaos. From a different perspective a similar conclusion was 
reached in [19] by Papenbrock and Weidenmtiller. As they state: "While the number of 
independent random variables decreases drastically as we follow this sequence, the com- 
plexity of the (fixed) matrices which support the random variables, increases even more. 
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In that sense, we can say that in the TBRE, chaos is largely due to the existence of (an 
incomplete set of) symmetries." 



6. Summary and future outlook 

In summary, we have introduced the embedded ensemble EGUE(2)-S'L'^(4) in this paper 
and our main emphasis has been in presenting analytical results. Our study is restricted to 
f/(f2) irreps of the type {4'",p}, p = 0, 1, 2 and 3. Using Eqs. (13)-(15) and the formulas 

for the functions P's, {H'^)"^''^'" , Q's and i?'s given in Tables 1 and 2, cross correlations 
in spectra with different (to, fm) irreps are studied with results presented in Fig. 2 and 
Section 5. See [2,9,5] for further discussion on the significance of cross correlations gen- 
erated by embedded ensembles (they will vanish for GE's). Elsewhere we will discuss the 
results for EGOE(2)-S'J7(4) and in the limit fl ^ oo the results for these two ensembles 
are expected to coincide except for a difference in scale factors. 

In future we also plan to investigate EGUE(2)-S'C/(4) for general U{Q) irreps for any 
TO and this is indeed feasible with the tabulations for sums of Racah coefficients given in 
[18]. Then it is possible to examine the extent to which EGUE(2)-S'C/(4), i.e. random 
interactions with S'C/(4) symmetry, carry the properties of Majorana or the C2 [S'C/(4)] op- 
erator. This study is being carried out and the results will be presented elsewhere. With 
this, it is possible to understand the role of random interactions in generating the differ- 
ences in the ground state structure of even-even and odd-odd N=Z nuclei. See [20] for 
a numerical random matrix study of N=Z nuclei. In addition, just as the pairing cor- 
relations in EGOE(l+2)-s have been investigated recently [6], it is possible to consider 
SU{n) D SU{'i), where SU{3) is Elliott's SU{3) algebra [21], and examine rotational 
collectivity with random interactions. To this end we plan to analyze in future expectation 
values of the quadratic Casimir invariant of SU (3) or equivalently that of quadrupole- 
quadrupole (Q.Q) operator over the EGOE(l+2)-S'C/(4) ensemble. Finally, going be- 
yond EGUE(2)-S'?7(4), it is both interesting and possible (by extending and applying the 
S'C/(4) D SUs{2) (g) SUt{2) Wigner-Racah algebra developed by Hecht and Pang [22]) 
to define and investigate, analytically, the ensemble with full SU{4) — ST symmetry. 
In principle, it is also possible to construct the to particle H matrix, which is SU (4) or 
SU{4:) — ST scalar, on a computer and analyze its properties numerically but this is for 
future. 
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